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Mirror dark matter provides a simple framework for which to explain the DAMA/Libra 
annual modulation signal consistently with the null results of the other direct detection 
experiments. The simplest possibility involves ordinary matter interacting with mirror 
dark matter via photon-mirror photon kinetic mixing of strength e ~ 10~^. We confirm 
that photon-mirror photon mixing of this magnitude is consistent with constraints from 
ordinary Big Bang nucleosynthesis as well as the more stringent constraints from cosmic 
microwave background measurements and large scale structure considerations. 
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A mirror sector of particles and forces can be well motivated from fundamental con- 
siderations in particle physics, since its existence allows for improper Lorentz symme- 
tries, such as space-time parity and time reversal, to be exact unbroken microscopic 
symmetries [1]. The idea is to introduce a hidden (mirror) sector of particles and forces, 
exactly duplicating the known particles and forces, except that in the mirror sector the 
roles of left and right chiral fields are interchanged. We shall denote the mirror particles 
with a prime ('). In such a theory, the mirror protons and nuclei are naturally dark, 
stable and massive, and provide an excellent candidate for dark matter consistent with 
all observations and experiments [2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16]. For a 
review, see e.g. ref. [17]. Dark matter from a generic hidden sector is also possible, see 
e.g. ref. [18] for a recent study. 

It has been shown in ref. [19], up-dating earlier studies[20], that the mirror dark mat- 
ter candidate is capable of explaining the positive dark matter signal obtained in the 
DAMA/Libra experiment [21], while also being consistent with the null results of the 
other direct detection experiments. The simplest possibility sees the mirror particles cou- 
pling to the ordinary particles via renormalizable photon-mirror photon kinetic mixing[22] 
(such mixing can also be induced radiatively if heavy particles exist charged under both 
ordinary and mirror C/(l)e^[23]): 

'-"mix 2 A"^ ^ ' 

where F^" = d^A' — d^A^ and F''^" = d^A'" — d'^A'^. This mixing enables mirror charged 
particles to couple to ordinary photons with charge eqe, where q — —1 for e', q — +1 for p' 
etc. The mirror dark matter interpretation of the DAMA/Libra experiment requires[19] 
e ~ 10~^, which is consistent with laboratory and astrophysical constraints [24]. 

The purpose of this note is to study the implications of such mixing for the early 
Universe. In particular, we will check that this kinetic mixing is consistent with constraints 
from ordinary Big Bang nucleosynthesis (BBN) as well as more stringent constraints from 
cosmic microwave background (CMB) and large scale structure (LSS) considerations. 

In the mirror dark matter scenario, it is assumed there is a temperature asymmetry 
(T' < T) between the ordinary and mirror radiation sectors in the early Universe due to 
some physics at early times (for specific models, see e.g. [25]). This is required in order to 
explain ordinary BBN, which suggests that T'/T ~ 0.6. In addition, several analyses[7, 8] 
based on numerical simulations of CMB and LSS suggest T'/T ~ 0.3. However, if photon- 
mirror photon kinetic mixing exists, it can potentially thermally populate the mirror 
sector. For example, Carlson and Clashow[26] derived the approximate bound of e ~ 
3 X 10~^ from requiring that the mirror sector does not come into thermal equilibrium 
with the ordinary sector, prior to BBN. The inferred value of e ~ 10~^ is consistent with 
this bound, so that we expect the kinetic mixing to populate the mirror sector, but with 
T' < T. Assuming an effective initial condition T' -C T, we can estimate the evolution of 
T'/T in the early Universe as a function of e, and thereby check the compatibility of the 
theory with the BBN and CMB/LSS constraints on T'/T. 

Photon-mirror photon kinetic mixing can populate the mirror sector in the early Uni- 
verse via the process e+e~ — > e'^e'~. This leads to the generation of energy density in the 
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mirror sector of: 



^ = rie+rie- {avM0iS) (2) 

where £ is the energy transferred in the process, Vm0i is the M0ller velocity (see e.g. 
ref.[27]), and n^- ~ ~ ^T^. 

It is useful to consider the quantity: p'/p, in order to cancel the time dependence 
due to the expansion of the Universe [recall p = Tr^f^T^/SO]. Using the time temperature 
relation: 

( = 0.39-/^^ (3) 

with g = 10.75 and Mpi ~ 1.22 x 10^^ MeV, we find that: 

d{p'/p) -ne-ne+{avM0i£) O.GMpi 



dT TT^gTySO v^TS 

Let us focus on {(tvm^i^)- This quantity is: 



(4) 



J l+e^i/T i+e£;2/r" Pl^ P2 

where we have neglected Pauli blocking effects. If one makes the simplifying assumption 
of using Maxwellian statistics instead of Fermi-Dirac statistics then one can show (see 
appendix) that in the massless electron limit: 

{avM0i£) = , (6) 



and Eq.(4) reduces to: 



d{p'/p) ^ (7) 



where 



Note that the e'"^ will thermalize with 7'. However, because most of the e'"^ are produced 
in the low T' ~ 5 MeV region, mirror weak interactions are too weak to significantly 
populate the i^'e^^^r [i-®- can easily verify a posteriori that the evolution of T'/T for the 

parameter space of interest is such that GpT'^ ^ o^Mpi ] • Thus to a good approximation 
the radiation content of the mirror sector consists of e'^, 7' leading to gf' = 11/2 and hence 
P'/P = {g'/gW^T^), with g'/g ^ 22/43. 



3 



Eq.(7) has the analytic solution: 



r ( g .V'^w 111/4 



(9) 



where we have assumed the initial condition T" = at T = Tj. 

Let us now include the effects of the electron mass. With non-zero electron mass, the 
evolution of T'/T cannot be solved analytically, but Eq.(4) can be solved numerically. 
Note that the number density is: 



'E^ - mlE 
+ exp{E/T) 

and, as we discuss in the appendix. 



1 roo JE^ - miE 



(11) 

where the cross section is: 

^ = f«Vl(. + 2m^r . (12) 

Numerically solving Eq.(4) with the above inputs (i.e. numerically solving the integrals 
Eq.(lO) and Eq.(ll) at each Temperature step), we find that^ 

e ~ 8.5 X 10-^° (^^y (13) 

where Xf is the final value (T — > 0) of x = T'/T. In figure 1, we plot the evolution of 
T'/T, for e = 8.5 X 10"^°. 

In deriving this result we have made several simplifying approximations. The most 
significant of these are the following: a) Using Maxwellian statistics instead of Fermi- 
Dirac statistics to simplify the estimate of {avM0i^)- Using Fermi-Dirac statistics should 
decrease the interaction rate by around 8% as discussed in the appendix, b) We have 
neglected Pauli blocking effects. Including Pauli blocking effects will slightly reduce the 
interaction rate since some of the e'^ states are filled thereby reducing the available phase 
space. We estimate that the effect of the reduction of the interaction rate due to Pauli 
blocking will be around ~ 10%. c) We have assumed that negligible t'e./i.r produced 
via mirror weak interactions from the e'^. Production of J^e.^x.r slightly decrease the 
T'/T ratio. The effect of this is equivalent to reducing the interaction rate by around ~ 
10%. Taking these effects into account, we revise Eq.(13) to: 

e = (1.0 ± 0.10) X 10"^ (5^)^ ■ (1^) 
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^For simplicity we have neglected the effect of heating of the photons via e+e~ annihilations. Note 
that the same effect occurs for the mirror photons which are heated by the annihilations of e'+e'", so 
that Xf is approximately unchanged by this effect. 
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Successful large scale structure studies [7, 8] suggest a rough bound on Xf of ~ 0.3. 
Our result, Eq.(14), then suggests a rough bound on e of e ~ 10^^. 
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Figure 1: Evolution of a; = T'/T for e = 8.5 x 10"^*^. The solid line is the numerical 
solution including the effects of the electron mass, while the dashed line is the analytic 
result [using Eq.(9], which holds in the massless electron limit. As expected the two 
solutions agree in the T ~ 1 MeV region, where the effects of the electron mass should 
be neghgible. 



In conclusion, previous work has shown that the mirror dark matter candidate can 
explain the DAMA/Libra annual modulation signal consistently with the null results of 
the other direct detection experiments provided that there exists photon-mirror photon 
kinetic mixing of strength e ~ 10~^. Here we have examined the implications of this 
kinetic mixing for early Universe cosmology, where we showed that it is consistent with 
constraints from ordinary BBN and CMB/LSS data. 



Appendix 



Here we shall examine the quantity {(jvm^i^) and derive Eq.(6) and Eq.(ll) used in 
our analysis. Following ref.[27], we have: 

, p. _ JavM0i{Ei + E2)e-^^/^e"^^/^d3pid=^p2 ,x 
{avM,it} - J^^^U^'^^^^^J^dF^;^, ^^^^ 
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where pi and p2 are the three-momenta and Ei and E2 the energies of the colhding 
particles in the cosmic comoving frame. Recall that E — Ei + E2 is the energy transfer 
per coUision. As elaborated in ref.[27], evaluation of these integrals can be facihtated by 
changing variables to E± = Ei± E2 and s — 2ml + 2E1E2 — 2piP2 cos 9. In terms of these 
variables the volume element becomes 

d^Pid^P2 = 277^ EiE2dE+dE_ds (16) 

and 

J avM0iEe-^'/^e-^''^d^pid'^p2 = 2t:'' j dE+E+ J dE^ J dsavM0iEiE2e-^+/^ 

(17) 

with integration region \E_\ < y^l — E\ — s, E^ > y^, s > Ami. Performing the 
E_ integration, we have: 

j avM0iSe-^'^''e-^'/^d'p,d'p2 = An'' J dsaF^l-^ J dE+e'^+Z^ El - s E+ 

(18) 

where aF — (7VM01E1E2 and F — \^ s{s — Ami). Also, as discussed in ref.[27] 

j e-^'/^e-'^^l^d^pid^p2 = [AnmlTK2{me/T)^ (19) 
where Ko is the modified Bessel function of order 2. Hence we see that 



(20) 

A 2 2 

In the me ^ limit, where a — ^3^^ , and using the dimensionless variable z = 
E+/y/s, we find: 

lavM0i£e-''^/^e-''^/^d'p,d'p2 = ^-^^^ dss^'^ dze-'^^'^ z^^^ 

= 128Q;2eV^r^/ (21) 

where 



Also, 

I e-^^/^e-^'/^d'p^d^P2 = [A7TmlTK2{me/T)f 

= 647r^r^ in the me ^ hmit. (23) 
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Thus we find: 

l^avM^iE) = -^r- ■ (24) 

Our results for {crvM0iS), Eq.(20) [or Eq.(24) for the rrie — > hmit], have assumed 
Maxwelhan distributions for the fermions to simphfy the integrals. In the rrie — ^ limit, 
it is possible to evaluate the integrals for the realistic case of Fermi-Dirac distributions. 
In which case, one finds: 

{avM0i£) = ^gy^ (25) 

where 



roo i-oo rV x'^ — z 

h — dz dx dy 

Jo JJz Jo 



X 



We find numerically that: 



^ Jo 1 + e^'+y 1 + e^'-y 

h = / dz dx dy y (26) 

Jo J-Ji Jo 1 + e^^y 1 + y 



h ^ 0.39, h ^ 0.84 =^ ^ ^ 0.46. (27) 



Thus, we see that the approximation of using Maxwellian statistics overestimates {(tvm^i^) 
by around 8%. 
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